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Abstract 

^-H I In Cator and Lopuhaa [3] an asymptotic expansion for the MCD estimators is established 

. in a very general framework. This expansion requires the existence and non-singularity of 

the derivative in a first-order Taylor expansion. In this paper, we prove the existence of this 
derivative for multivariate distributions that have a density and provide an explicit expression. 
■ Moreover, under suitable symmetry conditions on the density, we show that this derivative 

is non-singular. These symmetry conditions include the elliptically contoured multivariate 
location-scatter model, in which case we show that the minimum covariance determinant 
(MCD) estimators of multivariate location and covariance are asymptotically equivalent to 
a sum of independent identically distributed vector and matrix valued random elements, 
, respectively. This provides a proof of asymptotic normality and a precise description of the 

' limiting covariance structure for the MCD estimators. 

^ ■ 1 Introduction 

oo 

O 

I location and scatter parameters. These estimators, in particular the covariance estimator, also 

serve as robust plug-ins in other multivariate statistical techniques, such as principal component 
analysis O I18j , multivariate linear regression [H I17j , discriminant analysis [7] , factor analysis [13] , 
^ . canonical correlations [121 HO] and error- in- variables models [S] , among others (see also [B] for a 

I more extensive overview). For this reason, the distributional and the robustness properties of the 

MCD estimators are essential for conducting inference and perform robust estimation in several 
statistical models. 

The MCD estimators have the same high breakdown point as the minimum volume ellipsoid 
estimators (e.g., see [IllIT]). The asymptotic properties have first been studied by Butler, Davies 
and Jhun [2] in the framework of unimodal elliptically contoured densities, who showed that the 
MCD location estimator converges at ^/n-J:ate towards a normal distribution with mean equal 
to the MCD location functional. In the same framework, Croux and Haesbroeck [5] give the 
expression for the influence function of the MCD covariance functional and use this to compute 
limiting variances of the MCD covariance estimator. The asymptotic theory has been extended 
and generalized in Cator and Lopuhaa , who studied the MCD estimators and the corresponding 
functional in a very general framework. They establish an asymptotic expansion of the type 

1 " 

e„ -00 = -A'(0o)"'- V {"fix,, 9o) - E*(X„ Oo)) + or{n-^/^), (1.1) 
n ^ — ' 

i—l 

where On and Oq denote vectors consisting of the MCD estimators and the MCD functional at 
the underlying distribution, respectively, and ^(•,6'o) is a function that we will specify later on. 
In principle, from this expansion a central limit theorem for the MCD estimator can be derived. 



The MCD estimator [TFj is one of the most popular robust methods to estimate multivariate 
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However, the expansion requires the existence and non-singularity of A'{9q). Moreover, a more 
exphcit expression of its inverse is desirable from a practical point of view, since it determines the 
limiting variances. 

In this paper we show that A'(0o) exists as long as the underlying distribution P has a density 
/. Its expression given in Theorem 13.11 offers the possibility to estimate the limiting variances of 
the MCD estimators in any model where P has a density. We will also provide sufficient symmetry 
conditions on / for A'{9o) to be non-singular. This includes the special case of elliptically contoured 
densities 

fix) = det{Y.)-^/'^h{{x - fi)T.-\x - 

for which we show that the MCD location and covariance estimator are asymptotically equivalent 
to a sum of independent vector and matrix valued random elements, respectively. This exact 
expansion shows that at elliptically contoured densities the MCD location and MCD covariance 
estimator are asymptotically independent and yields an explicit central limit theorem for both 
MCD estimators separately, in such a way that the limiting covariances between elements of the 
location and covariance estimators can be obtained directly from the covariances between elements 
of the summands. Furthermore, the expansion for the MCD estimators is needed to obtain the 
limiting distribution of robustly re- weighted least squares estimators for (/i, S), if one uses the 
MCD estimators to assign the weights (see '12'). 

The paper is organized as follows. In Section [2] we define the MCD estimators and MCD 
functionals and discuss some results from ^ that are relevant for our setup. In Section [3] we 
establish the expression for the A'{9o) in terms of a linear mapping and show that this mapping 
is non-singular under suitable symmetry conditions. The special case of elliptically contoured 
densities is considered in Section [3J where we obtain an explicit expression of A'{9o)~^. From this 
we derive an asymptotic expansion for the estimators, prove asymptotic normality, and derive the 
influence function of the MCD functionals. As special cases we recover results from [2] and [4] 
under weaker conditions. 

All proofs have been postponed to an appendix at the end of the paper. 



2 Definition and preliminaries 

For a sample Xi,X2, . . . , Xn from a distribution P on M'' the MCD estimator is defined as follows. 
Fix a fraction < 7 < 1 and consider subsamples S C {Xi, . . . ,Xn} that contain /i„ > [717] 
points. Define a corresponding trimmed sample mean and sample covariance matrix by 

^■(*) = TT E 

(2^1) 

Cn{S) = I] " f^iSmx, - f^iS))'. 

Note that each subsample S determines an ellipsoid E{Tn{S),C„{S),r„{S)), where for each £ 
SJ', E symmetric positive definite, and p > 0, 

E{fi,j:,p) = {x e M'= : {x ^ nyT.-\x - fi) < p^}, (2.2) 

and 

?n{S) = inf {s > : P„ (-B(f„(5), Cn{S), s)) > 7} , (2.3) 

where P„ denotes the empirical measure corresponding to the sample. Let Sn be a subsample that 
minimizes det(C„(5)) over all subsamples of size hn > [^7], then the pair (X"^ (iSti) , Cyj (5*7^) ) is 
an MCD-estimator. Note that a minimizing subsample always exists, but it need not be unique. 
In [3] it is shown that a minimizing subsample Sn always has exactly [717] points and is contained 
in the ellipsoid E(Tn{Sn), Cn{Sn) Trn{Sn)) , which separates Sn from all other points in the sample. 
Note that in [2] (among others) one minimizes over subsamples of size \n^\ ■ This is somewhat 
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unnatural, since it may lead to subsamples S for which Pn{S) < 7. Moreover, it may lead to 
situations where the trimmed subsample does not contain the majority of the points, e.g., if 
7 = 1/2 and n is odd, then [njl = (n — l)/2. By considering subsamples S of size hn > ['1.7] in 
definition (|2.ip . we always have > 7 and for any 1/2 < 7 < 1, the subsample contains the 

majority of points. 

We define the MCD functionals in similar fashion. Define a trimmed mean and covariance as 
follows, 

^ \ \ (2.4) 
Cpi.^) - -TTlTi / ~ Tpm(x - Tp{<p)y^{x) P{dx). 



^ (j)dP 



and define 



rp{4,)=ini{s > ■.PiE{Tp{cl>),Cp{^),s)) > 7}. 

for measurable : R'^ ^ [0,1], such that j 4>dP > 7 and J \\x\\'^(p{x) P{dx) < 00. Note that for 
P — Pn and (j) = Is, we recover (|2.ip and (|2.3p . If <j)p minimizes dct(Cp(0)) over all considered 
above, then the pair {Tp{(j)p), Cp{(j)p)) is called an MCD functional. In [3] it is shown that such a 
(j)p always exists and a characterization of a minimizing </) is provided. From this characterization 
(Theorem 3.2 in ^3]) it follows that if P has a density, then 

(jyp^lEp and P{Ep)=-f, (2.5) 



where Ep = E{Tp{(t)p),Cp{<j)p),rp{<j)p)). This means that the MCD functional defined by ((2^ 
coincides with the definition through minimization over bounded Borel sets given in [5] . 

Throughout the paper we will assume that the MCD functional at P is uniquely defined, and 
we write (/xo,So) = {Tp{(l)p),Cp{(j)p)) and pq = rp((j)p). This holds, for instance, if P has a 
unimodal elliptically contoured density (see Theorem 1 in [^). We will also assume that P has a 
density / that satisfies the following condition: 

(B) / is continuous and strictly positive on a open neighborhood of the boundary of 
-B(^o,So,po)- 

In that case, it follows from Theorem 4.2 in [3] that 9n 6q with probability one, where 

On = (f„(5„), B„(5„),r„(5„)), with B„(5„)2 = C„(5„), 6) 
^'o (mo, To, Po) , with Tp = Eq- 

Moreover, Theorem 5.1 in |3] implies that expansion (|l.ip holds, where ^' = (^fi, ^'2, ^3), is 
defined as 

*l(2/, ^) = l{||G-i(j/-m)||<r}G"^(?/ - m) 

*2(j/, e) - lu\G-Hv-m)\\<r}{G-Hy - m){y - mYG-' ~ h) (2.7) 

^3(2/, 6*) = l{||G-i(a-m)ll<r} ~ 7, 

and A = (Ai, A2, A3), with 

A,(e) = j -^jiy, 0) P{dy), for 3 = 1, 2, 3, (2.8) 

for e = {m,G,r), with y,te'R'',r> 0, and G G PDS(fc). Here, PDS(fc) denotes the space of all 
positive definite symmetric k x k matrices. 
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3 Existence and non-singularity of A'{9q) 



Let ^0 — (mo: To, po) be the MCD functional at P. Due to the indicator function in the expression 
of ^{x,9o) it can be seen, that the existence of a derivative of A{9) at Oq cannot be expected in 
general if P does not satisfy condition (B). If P does satisfy (B), then the derivative will depend 
on the behavior of / on the boundary of E{^o, '^o, Po)- For p > and p G K'', define 

B{n,p) ^ {xeR'' : \\x- fi\\ < p} , 

and let aa denote the surface measure on the boundary dB{0,po). To describe the derivative of 
A{9) at 9o, we introduce the measure 

i^{duj) = det(ro)/(row + fio)ao{duj) for u e dB{0, po). (3.1) 

Here, ctq denotes the Lebesgue surface measure on 95(0, po)- 

Note that our parameter 6*0 is an element of 6 = M*^ x PDS(/c) x M. This means that the 
derivative of A at 9q, if it exists, can be described as a linear mapping on the tangent space of 
Q in 6*0, which is given by R*^ x S{k) x R. Here, S{k) denotes the space of all symmetric k x k 
matrices. The derivatives of Ai, A2 and A3 are given as linear mappings by the following theorem. 

Theorem 3.1. Suppose that P satisfies (B) and let the MCD Junctional 60 = (pq^Tq, po) be 
uniquely defined at P. For j = 1,2,3, the derivatives of Aj are given by the following linear 
mappings, with {h,A, s) € R'^ x S{k) x R; 

JdBo V Po ^PO / 

Kieo)ih,A,s) = -jir^'A + AT^')+ [ f^^yio!h^ -'i^o'A + AT,')u ^\ 

JdBo \ Po 2po / ^ ^ 
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A',(9o)(h,A,s)^ [^^+ ' ° ' °-^ + s]Hd^), 

JdBo V Po 2po 

where Bq = B{0, po) and ^{duj) is defined in (13. ip . 



Note that Theorem 13.11 also has practical implications. According to Theorem 5.1 in 3J, the 
MCD estimator 0„ = (p:„,r„,p„), represented as a vector, is asymptotically normal with mean 
zero and limiting variance given by the covariance matrix of Z = A'(^o)~^^(-^ij ^o)- This means 
that the expression for A'{9o) enables one to estimate the limiting variance of the MCD estimators 
in any model where P has a density, which goes far beyond the traditional elliptically contoured 
densities. An estimate for A'{0q) can obtained by plugging in the estimate 0„ for Oq, replacing ctq 
by the surface measure (t„ on dB{0, pn), and using a nonparametric estimate for the density / on 
the boundary of -6(0, pn), e.g., histogram or kernel type estimates. As long as the estimate A'(6'„) 
turns out to be non-singular, the limiting covariance matrix of \fn[()n — ^0) can be estimated by 
the sample covariance of the Zi = A'(^?„)^^5'(A'i, 0„). 



We proceed by finding sufficient conditions for A' (6*0) to be non-singular. We would have non- 
singularity, if for all (/i. A, s) e R'= x S(k) x R, A'{9o){h, A,s) = implies that {h, A, s) = (0, 0, 0). 
From the expressions in Theorem 13.11 it can be seen that this cannot be expected without further 
assumptions on /. Suitable symmetry assumptions on / will simplify the expressions for the 
derivative, in which case non-singularity can be established. First note that without any further 
assumptions on /, we do always have the following result. 

Lemma 3.1. Suppose that P satisfies (B) and let the MCD functional 6a = (/iq, Fq, po) be uniquely 
defined at P. Let A be defined by (12. 8p and suppose A'{0o){h,A,s) = 0. Then Tt{Tq^A) = 0. 
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Next, we consider the case where the density / is point symmetric with respect to the center 
of £;(/xo,So,po), i-e., 

f{-ToCJ + fio) = f{To^ + f^o), foTCjedBiO,po). (3.2) 
In that case we have the foUowing lemma. 

Lemma 3.2. Suppose P satisfies (B) and (13. 2p . and let the MCD functional 6q = (fio,To, po) be 
uniquely defined at P. Let A be defined by (j2.8p and suppose that A'{9o){h, A, s) — 0. Then 
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uj'{Tf^'A + ATQ')ujiy{duj), (3.3) 



2^01^0 

where Bq = -6(0, po) o-nd vq = vIOBq), with v defined in (13. ip . If, in addition, for all i = 1, 2 . . . , fc 

uj1v{duj) ^-IPo, (3.4) 

dBa 

then h — 0. 

Point symmetry will not be sufficient to conclude that A = from A' {dQ){h, A, s) = 0. The 
slightly stronger condition of half-space symmetry will suffice, i.e., 

f{TQUj(_,) + fio) = /(row + /io), where = (wi, . . . , -w^, cj^+i, . . . , w/c), (3.5) 

for all i — 1,2, ... ,k and uj G dB{0, po). To describe sufficient conditions for non-singularity, we 
define the matrix Af with elements 

Mij ^ 1 ijj'luj'j v{du) - — [ uj1v{du) I w| i/(du;) - 27pol{i=j}, (3.6) 

JdB{0,po) ^0 JdB{0,po) JdB{0,po) 

for i,j = 1,2, ... ,k, where i^a = v{dB{Q, po)) and v{du)) is defined by p.ip . We then have the 
following theorem. 

Theorem 3.2. Suppose that P satisfies (B) and (13. 5p . and let the MCD functional 9q = {p,Q, Tq, po) 
be uniquely defined at P. Suppose that (|3.4p holds, that for all i,j = 1,2 ... ,k with i ^ j, 

Lufui'^ v{dw) ^ 7po, (3.7) 

'9-8(0, po) 

where v is defined in (13. ip . and that the matrix M defined in (j3.6p is such that for any x € M.^ , 

Ma; = and xi ^ h Xfc = a; = 0. (3.8) 

Then, for do = [fiQ, Tq, po), the derivative A' (do) is non-singular as a linear map on R'^ x S{k) x M. 

Example of densities that satisfy p.Sp are elliptically contoured densities. However, also aSine 
transformations of densities that have independent marginal densities that are symmetric around 
zero, i.e., 

f{x) = g{T^^{x - p)), where ^(xi, ...,Xk)^ 9i{xi) ■ ■ ■ gk{xk) and gi(xi) = gi{-Xi), 
satisfy (13. 5p . 



5 



4 Elliptically contoured densities 

Suppose that P has an ehiptically contoured density, i.e., 

f{x) =^det{T.)-^/'^h{{x- ^i)J:-\x ~ H)) where A* e M*-',!] e PDS(/fc), (4.1) 

and h : [0, oo) [0, oo) is decreasing so that P is unimodal. In this case, it follows from the 
characterization for the </> function that minimizes det{Cp{(j))) (see Theorem 3.2 in [3]), that our 
definition of the MCD functional coincides with the one in [2] , who show that the MCD functionals 
are unique: 

Mo=/^, So=«(7)'S, and pl^^dl^, (4.2) 



where 

and where r(7) is determined by 

h{r'^)r^-^ dr = -f. (4.4) 



r(fc/2) 7o 



The next proposition shows that for elliptically contoured densities the derivative A' (6*0) exists 
and is non-singular. 

Proposition 4.1. Let P have an elliptically contoured density as defined in (|4.ip with h non- 
increasing such that P is unimodal. Then all conditions of Theorem \3.2\ are satisfied. 



We proceed by obtaining asymptotic expansions for the MCD estimators in the case of ellipti- 
cally contoured densities. Because the estimators are afhne equivariant, it suffices to consider the 
spherically symmetric case (^, S) = (0, J). The next theorem provides the expressions for A'(6'o) 
and its inverse at spherically symmetric densities. 

Theorem 4.1. Let P have an spherically symmetric density f{x) — h{\\x\\'^) with h decreasing 
such that P is unimodal. Let r = r{^) and a = a(7) he defined in (|4.4p and (|4.3p . respectively, 
and let D — K'[9q), for Oq — {fio,TQ, po) = (0, a/, r/a). Then the linear mapping D is given by 



D2 



{h,A,s) i-> Pih 

(/i. A, s) ^ P2A + (5zTr{A) ■ I + ff^s ■ I 
(h, A, s) ^ /?5Tr(A) + Pes, 



and the inverse linear mapping D™^ is given by 
[D'^^]^:{g,B,t)^P^'g 

[D ]^:{g,B,t)^P, B + ^^-^^-^Tm-I+^^t-I 
[i^-]3:(^^,5,t)-|^Tr(i?)-^^%^^ 



where 



/3i-^(f^o-7) <0, /34 = ^^o-^o, 
ak{k + 2) a ' fca ' 
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with Bq — i?(0, po) fl^fl 



An immediate consequence of Theorem 14.11 is the next corollary, which shows that the MCD 
estimators of location and covariance are asymptotically equivalent to a sum of independent iden- 
tically distributed vector and matrix valued random elements, respectively. 

Corollary 4.1. Suppose that P has a spherically symmetric density f{x) — with h de- 

creasing such that P is unimodal. Let r — r(7) and a — 01(7) he defined in (|4.4p and (|4.3|) . 
respectively. Then for n ^ 00, 



1 " 

y/np-n = —i='^'n:t{\\x,\\<r}Xi + op(l); 



'n . 



1 " 

V^(I]„ - a^I) = ^ XI [^{\\X^\\<r} ('^1 • ^ + ^^2\\X^\\^ ■ I + K3^.X;) + K4 • /] + Op(l); 
7" 1 " 

(Pn - -) = ^ X! [•^ll{l|X.|l<'-}ll^ill^ + •^2l{||X,||<r} + A3] + Op(l), 

2 — 1 



where tt = — (a/3i) ^ and 



al32 + 27 2 - fca^ 



f^i — — — ; ^ — , '^a — '«4 — , 

K7 Kjap2 ap2 K 

" "2fc^' " 2fc^ " ^' " ^ ^ 2^3 



wif/i /3i , f32 and f3e defined in Theorem \4.1 



We proceed by obtaining the limit distribution of the MCD estimators. To describe the limiting 
distribution of a random matrix we use the operator vec(-) which stacks the columns of a matrix 
M on top of each other, i.e., 

vec(M) = {Mil, Mik, . . . , Mfci, . . . , Mkk)' ■ 

We will also need the commutation matrix Ck,k, which is a fc^ x fc^ matrix consisting of fc x fc 
blocks: Cfe.fc = (Aij)^^^]^, where each (i, j)-th block is equal to a fc x fc-matrix Aji, which is 1 at 
entry (j, i) and everywhere else. Finally, for matrices M and N, the Kronecker product M (g) N 
is a fc^ x fc^ matrix consisting of fc x fc blocks, with the (i, j)-th block equal to niijN. 

Theorem 4.2. Suppose that P has a spherically symmetric density f(x) = ft.(||a;|p) with h de- 
creasing such that P is unimodal. Let r = r(7) and a = a(7) be defined in (j4.4p and (j4.3p . 
respectively. Let fin, S„ and pn the MCD estimators. Then 

(i) fin and (E„,p„) are asymptotically independent, the diagonal elements oj Yin asymptot- 
ically independent from the off-diagonal elements and pn, and the off-diagonal elements of 
Tin are asymptotically mutually independent; 

(a) ^/npn is asymptotically normal with mean zero and covariance matrix tI, where 

{kja — rvo)"^ ' 

where i>o is defined in Theorem\4-.l\ 
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(Hi) Y^(vec(S„) — Q!^vec(/)) is asymptotically normal with mean zero and covariance matrix 

<Ji{I + Ck,k){I (g>I) + cr2vec(/)vec(/)', 

where 

2 1 .,4 jr'^ ~2kjr'^a'^ + k'^ja'^ + 2kr'^a'^ - r'^ 

'^^ = -k^^ + ]^^M\\X.\\<r}\\X,\\ — 

where K3 is defined in Corollary \4-3[ 
(iv) ^yn{pn — r I a) is asymptotically normal with mean zero and variance 

2 2,1? ii4 2fc^fora^7 — I'Qfcr^Q!^ + 4fc^7^Q!® — 4fcr'^z^oQ!^7 + ?'^^'o 

(7 = AiEl{||Xi||<r} ^1 H ., 2 6 2 ' 

where Ai is defined in Corollary \4-.l\ and vq is defined in Theorem \4.1\ 

With Theorem I4.2f i) we recover Theorem 4 m [2j. Note however, that the assumption of h 
being difFerentiable (see [2]) is not required in our approach. Furthermore, it can be seen from the 
expression of the hmiting variance of E„ that in the spherically symmetric case: 

^/nY^njj N{0,ai) 

^iS..-.^ j^^lU (7. 2.l+(7jj' 

for i,j l,2,...,fc._ 

Because /i„ and S„ are affine equivariant the limiting distributions for the MCD estimators in 
the case of general e M'^ and S € PDS(/c) can be obtained easily. Suppose that Xi, . . . , Xn are 
independent with density 

fix) = Acm-^'^h{{x - fifj:-\x - At)). 

Because of affine equivariance it follows immediately that ^yn{p,n—fi) is asymptotically normal with 
zero mean and covariance matrix T{tI)T ~ tS, where = E. Similarly ^/ri{vec(T}n) — a^vec(S)) 
is asymptotically normal with mean zero and covariance matrix Evec(FMF)vec(FMF)', where M 
is the random matrix with Evec(M)vec(M)' = ai{I + Ck,k) + (72vec(/)vec)/)'. It follows from 
Lemma 5.2 in [9J, that 

Evec(FMF)vec(FAfF)' = (7i(/ + Ck,k){^ «) S) + (72vec(S)vec(S)'. 

This means that we have the following general corollary of Theorem 14.21 

Corollary 4.2. Suppose that Xi, . . . , Xn are independent with an elliptical contoured density 

f{x) = dei{i:)-^l^h{{x ~ ^i)'Y,-^l'\x - ^)) , ^ e R^ S e PDS(A:), 

where h : [0,oo) [0,oo) is non-increasing such that f is unimodal. Let (/i„,I]„) be the MCD- 
estimators. Then 'fin and E„ are asymptotically independent, y/n(jln — /i) has a limiting normal 
distribution with zero mean and covariance matrix tYi and •\/n(vec(S„) — (af'^vec(S)) has a limiting 
normal distribution with zero mean and covariance moirix (7i(/ + CA;^fe)(E(8) S) + (72vec(E)vec(I])', 
where t, ai and oi are given in Theorem\ 
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Another corollary of Theorem 14.11 is the expression for the influence function of the MCD 
functional. The influence function of a functional &{■) at P is defined as 

IF(., e, P) = hnr ^ii^~e)P + eS.)-eiP) 

sio e 

if this limit exists, where 5x is the Dirac measure at a; € M*^. Denote by /i(-P) = Tp{(t)p), 
S(P) = Cp{(j>p), and p{P) — rp{(t)p) the MCD functionals at distribution P. We then have the 
following corollary. 

Corollary 4.3. Suppose that P has a spherically symmetric density f{x) — with h de- 

creasing such that P is unimodal. Let r = r(7) and a — a^-y) be defined in (j4.4p and (j4.3p . 
respectively. Then, for a; G M*^ such that \\x\\ ^ r, the influence functions of the functionals IJ.{P), 
S(P) and p{P) are given by 

IF(a;;^, P) = -kI {\\x\\<r}X 

IF(x; S, P) ~ l{||^||<r} (ki • / + K2||2:|P • / + ksX-x') + K4 • / 
IF(a;; p, P) = Ail{||^||<r}||a:||^ + \2'^{\\x\\<r} + A3 

where tt, Ki, K2, K3, K4 and Ai, A2, A3 are defined in Corollary \4-.l\ 

Clearly, all the expressions in Corollary 14. 31 are bounded uniformly for ||a;|| ^ '''(7). For a; S M'^ 
with ||a;|| ~ r{'-f), it is not clear whether the limit in (|4.5p exists, not even in the case of a unimodal 
spherically symmetric density. As a special case of Corollary 14.31 we recover Theorem 1 in [4] . 
However, we do not need the assumption that h is differentiable (see [3]). In order to see that our 
expressions coincide with the ones in [3], note that their quantities g, a, g^, and c^, correspond to 
our h, 1 —7, r(7)^, and l/a(7)^, respectively, and that they consider the Fisher consistent version 
of the covariance functional, i.e., Cq x S(P). Moreover, their expression bi — kb2 is simply equal 
to 1. For further discussion on IF(x; S, P) at spherically symmetric densities and corresponding 
graphs, we refer to [4]. 



5 Appendix 

5.1 Proofs for Section [3] 

The proof of Theorem 13.11 requires the following lemma, which helps to describe the derivative 
of A when (/io, Fq, s) — (0, /, r), in terms of a linear mapping. Let M{k) be the space of all fc x fc 
matrices. 

Lemma 5.1. Let r > and cj) : H.'^' — > M™, which is continuous on dB{0,r). Define the mapping 
L:R^ X M (fc) X M ^ M™ by 

L{h,A,s) = / (t)iy)dy. 

JE{h,{I+A){I+A)\r+s) 

Then, the derivative of L at (ho, Aq, sq) — (0,0,0), is given by the continuous linear mapping 
L'{0,0,0){h,A,s)^ [ (^ + ^^:iA±^ + ,)^^u;)ao{dio), 

with {h, A, s) e M*^ X M (fc) x M. 

Proof: The derivative can be found as the sum of the derivatives of 

/ '/'(y)dy, / '/'(y)dy, and / (j){y)dy. (5.1) 

JB(h,r) J E{0,{I+A)(I+A)' ,r) J B(0,r+s) 
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For the first integral, consider 

/ '^(y)dy- / 0(2/) dy= / (lij(ft.r) - li3(o,r)) 0(2/)d2/, 

JB(h,r) J B(0,r) J 

for 0. In first order this reduces to integration over dB{0,r). Let uj £ dB{Q,r), let 

V = {1 + S)uj E dB{h,r), and let a denote the angle between ui and h. Then the law of cosines 
yields that 

= \\vf + \\hf - 2\\v\\ ■ \\h\\ cosa = (^\\v\\ - jj^^ + \\h\f{sinaf. 

Since = r, in first order we find = (1 + 6)r'^ — oj'h, or S = (cj'/i)/r^. This means that for 
each u e dB{0,r), the length over which we integrate (j){uj) is — ||a;|| = S\\u}\\ = {ui'h)/r. Since 
(j> is continuous at u; e dB{0, r), we get, for 0, 

^{y)dy^ ( (j){y)dy^ f — ^{u) ao{duj) + o{\\h\\). 

B{h,r) J B{<d,r) JdB(0,r) ^ 

For the second integral in (|5.ip we consider 

(l£(0,(7+A)(/+A)',r) - Is(O.r)) (f>{y)<iy, 

for ^ 0, which reduces to integration over uj e dB{0,r). Let w = (1 + S)lu be such that 
||(/ + = r. Then 

2 

Since, for ^ 0, we have (/ + A')-\l + A)-^ = I - A - A' + 0{\\A\\^), it follows that 

This means that for each ui £ dB{0, r) the length over which we integrate (f>{i^) is 

I, uj'(A + A')uj „ ,,, , ,1 J, 

Ml - ll^ll = SM - + 0{\\Af). 

This implies that 

(lB(0,(/+A)(/+A)',r) - li3(0,r)) 0(2/)d2/ = / ^ '''^ ^ '/'(^) ^0 (t^^) + o( H ) ■ 



dB{a,r) 2r 



Finally, for the third integral in (|5.1I) we obtain 



(ls(0,r+s) - lB(0,r)) 0(y) = S / f7o(rfw) + o(s). 

Summing the three linear mappings, yields the desired result. ■ 

Proof of Theorem I3.lt First note that everything can be rescaled to the situation with iiq = 
and Fq = /, i.e., for any function g{y), we have 

'iE(po+h.{ro+A)iro+Ay.po+s){y)9(y)P{dy) 

r (5.2) 

det(Fo) / lj^0i^^j^^y^p^^^^{z)g{TQZ + ^ia)f{Toz + ^io)dz, 
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where h^T^^h and A = Tq ^A. To compute A'^{0o), take g{y) 1 in (|5.2p and for r/ = (/i, vl, ; 
(0, 0, 0), consider 



^siOo+v) - A3(6'o) = y T^E{f,o+h,{ro+A){ro+Ay,po+s){y)P{dy) - y l£(^o.So,po)(y) -^('^2/) 

= det(ro) J (l£;(X,(/+l)(/+l)',p„+s)W-l£(oj.Po)(^)) /(ro2 + Mo)d2 
= i'(0,0,0)(/i,I, s) + o(||(/i,A s)||), 
by taking (^(z) = det(ro)/(ro2; + /io) in Lemma [5TT] We conclude that 

A!,(0o) = det(ro) / ( + -'i^o'^ + A'T-')^ ^ \ ^^^^^ ^ 
JdB„ \ Po 2po J 

For the location functional, with 6 — (m, G, r), we have 

^ G-\y-m)Pidy)] + I Y,\y-i^,)P{dv)\ 

/ „ „ C''' \JE(m,G^,r) J 



-B(aio,So,po) / 



(5.3) 

The first term on the right hand side of (|5.3p can be decomposed as 



d_ 
00 



Eo 



d_ 

86 



T^\y-m)P{dy) 

Eq 



(5.4) 



where Eq = £'(/io, So, Po)- Because of (|2.4p and (|2.5p . it follows that 



P(£;(/io,So,Po)) =7 and / (y - ^lo) P (dy) = , (5.5) 

-'£;(/Jo,So,po) 



so that the first derivative in (|5.4p is equal to zero. To determine the second derivative in (|5.4p 
write 

/ To i(y - A*o - /i) P{dy) - f T^\y - a*o) P{dy) = -7ro-'/i, 

J_B(po,So,po) J_E(po,Eo,po) 



which yields 



To^^y - m) P(dy) 

-E(po,So,po) y 



= -7ro '/i. 

= 00 

For the second term on the right hand of (|5.3p . for (ft,, A, s) — * (0, 0, 0), consider 

J '^E{,,o+h,{ro+A}{ro+Ay.po+s}{y)'^o^{y - Mo) -P(rf2/) - J '^E{,,o,^o,Po)iy)'^o^iy - Mo) P{dy) 

= det(ro) y" (lB(7i,(/+A)(/+A)',po+s)(^) - l£(o,/.po)(2)) ^/(^o^ + Mo) dz 
= L'(Q,0,0)(h, A,s)+o{\\{h,A,s)\\), 
by taking (/)(z) = det(ro)2:/(roz + /ig) in Lemma [5. II It follows that 



d_ 
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T,\y ~ ^io) Pidy)] 

E(m,G^,r) J 

det(ro) / f ^ + + + s] com. + Mo)ao(rf-) 

JdB„ \ Po 2po y 

det(ro) / °— + ^ ° ^ ^ + s c^/(roco + /^o)^o(d^)- 
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We conclude that 

A'i(^o) = -7r(7'/i + det(ro) 
Then, similar to (I5.3I1. we have 



dBa \ PO 



2po 



s ujfiTouj + /io)o-o(dw). 



-B(aio,So,po) 



+ 
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E{m.,G^,r) 



G~\y^m){y-mYG-' -l\ P{dy) 



(5.6) 



The first term in (|5.6p can be decomposed as 

G-i(2/-/io)(2/-/^o)'ro-'-/] P{dy) 

T^\y - m){y - ^i^yV^^ - l] P{dy) 



d_ 
89 



d_ 

m 



d_ 



+ 



d_ 

de 



E{na.T,a,pa) 



■E(mo,So:Po) 



(5.7) 



T^\V - fio){y - m)'ro-i - /J P{dy)j 
To\y - Mo)(2/ - Mo)'G-i _ ^] p(^y) 

For the first term in (|5.7p . for ||^|| 0, consider 

(Fo + A)-'iy - fio){y - ^ioyTo' - l] P{dy) 

T-\y - /io)(2/ - Mo)'r„ ' - /] P{dy) 
{I - T^'A)T^\y ~ ^o)iy - Mo)'ro ' - 1] Hdy) + oi\\A\\) 

r^\y - ^^o){y - Mo)'ro ' P{dy) + o{\\A\\) = -^T^'A + o{\\A\\), 



E(no.T,o,po) 



' E{po,^o,Po) 



'■E(po,So:Po) 

where in the last two steps we use 

Ta\y - fio){y - Mo)'ro"' P{dy) = jl, 

I E(tiQ,Y,Q,po) 

which follows from (|2.4p . For the second term in (|5.7p . consider 



(5.8) 



-E(po,So,po) 



r^\y-pio-h)(y-fioYro'-I P{dy) 



B(po,So,po) 



ro'(y-/^o)(y-Mo)'ro'-/ P{dy) 



{y ^ fioYro' P{dy) ^ 0, 



-E(po,So,po) 

where we use (|5.5p and (|5.8p . Because G and Fq are symmetric, the last two terms in (|5.7p are the 
transpose of the first two terms in (|5.7p . This leads to the following derivative for the first term 
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in (ESI): 
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G-^{y-m){y-m)'G- 



P{dy) 



' E{fj.o,'So,po) 

For the second term on the right hand of (|5.6p . for (h, A, s) — > (0, 0, 0), consider 

^E{^,o+h,{^o+A){^o+Ay.po+s){y) [^o^iv ^ f^o){y - f^oY^o^ - l] P{dy) 

- J l£(po.Eo.po)(2/) [^o\y - Mo)(y - Mo)'ro ' - /] P{dy) 

= det(ro) J (l£;(X,(/+l)(/+l)',po+s)(^)-ls(o./,Po)(^)) [zz' - I]f{Toz + fio)dz 
= L'{0,0,0)(h,A,s) + o{\\{h,A,s)\\), 
by taking — det(ro)[zz' — I]f{TQZ + /io) in Lemma [01 It foUows that 



d_ 
80 



ro\2/-Aio)(y-/^o)'ro^ -/ P{dy) 



det(ro)/ 

JdBo \ PO 



Lo'h uj'{A + A')uj 



2po 



det(ro) 



We conclude that 



dBo 



TqI/i lo'{T-'^A + A'T-^)ijj 



Po 



2po 



+ s I {ujuj' - I) /{Tquj + fio)aQ{duj). 



A',{0o)ih,A,s) = -~jiT^'A+A'T,')+ f (- 

JdBo \ 



PO 

Noting that we take A symmetric, this finishes the proof. 
Proof of Lemma I3.lt From Theorem 13. II 
= A'M{h,A,s) 



2po 



= -liV^^A + AFo-i) + / °— + ^ ° ^ + s iuju' - I) v{duo), 



where Bo = i3(0, po)- Taking traces yields 



Jaso V Po ^Po / 

Because A3(6'o) = 0, if follows from Theorem 13.11 that the second term on the right hand side is 
zero, which proves the lemma. ■ 



Proof of Lemma [HI} If / satisfies then 

LOiUOjLOjn v{duo) ~ uJi v{dLo) — 0, for all i, j, m = 1, 2, . 

) JdBo 

where v{du}) is defined by p.ip . Hence, from Theorem 13. II we get 

= K'^{eo) I i^'{T-^A + AT-^)uj u{dLo) + si'idBo), 

2po JaBo 



, k, 



(5.9) 
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which yields the first statement. Moreover, (|5.9p and Theorem 13.11 also yield that 

O = A[{0o) = -ir-^h+ f ^° ^ ujv{duj) = — ( [ ujLu'i^{duj)^-jpal]ra^h. 

JdBa PO PO \JdBa / 

Because / satisfies p.2p , the matrix on the right hand side is a diagonal matrix with elements 

uj'^ iy{duj) - 7po- 



Therefore, from ([23) it follows that h = 0. 

The proof of Theorem 13.21 requires the following lemma. 



Lemma 5.2. Let A be a k x k matrix and suppose that TA + AT = 0, for some k x k positive 
definite symmetric matrix T. Then A — 0. 

Proof: Since F is positive definite symmetric, there exists a basis of eigenvectors of F. Choose 
V an eigenvector with eigenvalue A > 0. Then r{Av) + X{Av) = 0. This means that either Av 
is an eigenvector of F with eigenvalue — A < 0, which is impossible since F is positive definite, or 
Av = 0. This holds for all eigenvectors of F, and therefore A = Q. ■ 

Proof of Theorem 13. 2t Suppose that A'j{9Q){h, A, s) — for j = 1,2,3, then it suffices to 
show that {h,A,s) = (0,0,0). Because (13.51) implies p.2p . it follows from Lemma [3.21 that h — 0. 
Furthermore, A2(0o) = and A'^{Oo) = imply that 

O = A'2(0o) = ~75+ / (^^ + s]ujuj'j^{duj), (5.10) 
JdBo V 2po / 

where Bq = B{0,po) and S — Fq + ATq^ is symmetric. Condition (|3.5p implies that 

^j^^m '^{doj), for m = n;i=j, 

I 

^mt^« i'(rfw), for mj^n; {i,j} = {m, n} 



UJiUJjUJmUJn v{dljj) = < 

dBo 



dBo 
dBo 

0, otherwise. 



where v{du)) is defined by (|3.ip . Consider the (m, n)-th element of equation (|5.10p for m ^ n. 
Then it follows from (|5.1ip that 

= -27Po<S'„i„ + 25,„„ / LOm^n ^{duj) ^"^i I ^m^n v{duj) - 7po I '5'rrm- 
JdBo \JdBo J 

The factor in front of Smn is non-zero by assumption (|3.7p . so that Smn = for all m ^ n. Finally, 
consider the (m, m)-th element of (|5.10p and insert p.3p . which is obtained from Lemma |3^ Then 
we get 

2, ,2 ' ■ ^ ^ 



= -27Po<S'm„i + / ' 

\JdBo 



\ ujf v{duj) j 




JdBo J 


dBo / 



The right hand side is of the form Mx, where x = diag(S') and M is defined in (j3.6p . However, 
since Tr(S') = according to Lemma l3Tl from (|3.8p we conclude Smm = for all to = 1, 2, . . . , A:. 
It follows that S — T'^^A + ^Fq"^ = 0, and consequently, by p.3p . we have s = 0. Furthermore, 
from Lemma [5.21 we conclude that A — 0. ■ 
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5.2 Proofs for Section [3] 

Proof of Proposition [47ll Conditions (B) and (|3.5p arc immediate. From (|4.2p . wc find that 
/(Tquj + /io) is constant on dBo: 

fiToUj + Aio) = det(E)-i/2/,(^(^)2|,^|,2) ^ det(ro)-i/i(r(7)2), 

and 

lyidio) = a{j)''h{r{j)^)ao{dLo), (5.12) 
for u! € 9Bo- One can easily check that for all i, j = 1,2, . . . , k (e.g., see Lemma 1 in |10) ) 

1 r 9-^^/"^ 
-" + ^^«' Mtr2)r(./2) '"'''-"'"''^'''''°' 



Because h is decreasing and non-constant on [0,r(7)), conditions (|3.4p and (|3.7p are fulfilled: 

o fe/2 /•i-(7) ^tt'^/S 



Finally, from the equations above, it follows that the matrix M defined in (|3.6p can be decomposed 
as M = cil + C2II', where 1 = (1, 1, ... , 1)', and 

= fe(fc + 2)r(fc/2) '('^(")')^(^^'^" - 
- - fc^(fc + 2)r(fc/2) "(^(^^')^(^^ '^^ 

Because Mx = cix + C2(a;i + • • • + 2:^)1, it follows that, if + • • • + Xfe = 0, Mx = implies x — 
as long as ci 7^ 0, i.e., 

fc(fer2)r(./2) ^(-(^)^)-(^)''"^^^-(^)^' 

which follows from (|5.14p . ■ 
Proof of Theorem HU From (|5T^ . it follows that 

LUiLOj ^{duj) = 0, for i ^ j 

(5.15) 



uoi v{duj) — and / ui.iOjjUJk ^{diu) = 0, for all k, 



I dBo 

Hence, from Theorem 13. II we find 



A[{9o) : {h,A,s) ^ -jr^^h+ — [ cj^Tq^/ii^ 

PO J dBo 
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where according to (|5.13p and (|5.14p . 

kapo Jqbo a a \ k I 

Next, consider A3(6'o). From (|5.15p we find 

^3(^0) = — oj'r^^Aujiyiduj) + si^Q. 

PO JdBo 

From (|5.13p . the first term on the right hand side is 

— / uj'Auv{duj) = -^ [ \\uofTT{A)v{duo)^f^TT{A). 
This means that A!j(6'o) = /35Tr(A) + /^gs. Finally, from (jl^D and dEUD, 

Ai(^o) = -7(r„-M + Ar„-i) + / H i^o' A + AT,')oj ^ \ / , _ ^ ^^^^ 

JdBo V ^PO / ^ ^ 

= —A — AU9o)-I-\ / (uj' Auj)u!uj' iy{duj) + s [ uuj' u{dhj) 

" aPO JdBo JdBo 

= -'^A - [P^Tv{A) + ^gs) • / + — / {JAuo)loJ v{dLo) + ■ si. 
a apo JdBo k 

Consider the (m, n)-th element of the third integral on the right hand side. From (|5.11l) and (|5.13p . 
it follows that this integral is equal to 

1 ^ ^ f 

Aij / UOiUJjUJrnt^ri v{duj) = , f"^" {Tt {A)t Un=n} + 2A„„) , 

"^0 7^1 J^l JdBo "^(fc + 2) ^ 

which means that 

— / {JAu^)u:u' v{du:) = (Tr(A) ■I + 2A). 

"PO JdBo "«^(«^ + 2) 

Summarizing, in the expression of A2(0o), the coefficient of A is 

^ 2pluo _ 27 
'^^^ a/c(fc + 2) a ' 

the coefficient of Tr(^) • / is 

^ pIvq _ Poi^ 
afc(A; + 2) ka ' 

and the coefficient of si is 

/, Po 

Pi = -^^Q - t^O- 

From (|5.13p and (|5.14p . it can be seen that /32 < 0. 

To determine the expression of the inverse mapping, put D{h,A,s) — {g,B,t) and solve for 
{h,A,s). For the vector valued component of D, we have g = Di{h,A,s) = (3ih. Since /3i < 0, 
this immediately gives h = Pi^g. For the remaining mappings put 

B = D^ih, A, s) - 132A + I33Tt{A) -1 + 1348 -I 

t^D3{h,A,s)^l35MA)+l36S. ^' ' 
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By taking traces in the first equation we can solve for Tr(^) and s: 

cTr(A) = peMB) - kp^t 

cs = i(32 + k(33)t-f35TT{B), 

wlrere c = (32(36 + kPsPe - fc/34/35 = -2-^(3^^/ a. Since /32 < and /Je > 0, from (|57T6)) and (|5l7ll it 
follows that 

A = (3^^ {B - /33Tr(^) ■ I - p^s ■ I) 

- (3^^B ~ ^ (PeTiiB) - k(3it) -I-^ (-/35Tr(B) + {(32 + fc/Ss)^ • / 

CP2 C(32 

and 

, ^ -^Tr(i3) + h±]lht = ^Tr(i?) - ^(^l±^t. ■ 
c c 27/36 27/36 

Proof of Corollary 14. It Since A(6'o)~^ is a linear mapping and E^'(Xi,6'o) = 0, we obtain 
from (frT|l . 



n 



- V A'(0o)"'^(^., ^o) + ov{n-^l^), (5.18) 
n ^ — ' 

where 'I' is defined in (12. 7p . In particular, we have 

1 " 

1 " 



(5.19) 



/n 

i— 1 



According to (|4.2p . = (a'DiToiPo) = (0, a/, r/a), so that 

*i(a;,6'o) = l{||2,||<.r}a"^x, 

*2(a;, 0o) = l{|b||<r} (a^'xx' - /) , (5.20) 
*3(a;,6lo) = l{||x||<r} -7- 

Insert go ^'i(a;, ^o), -Sq = 5*2(2;, 6*0), and to = ^3(2^1 ^0) in the expressions for £''"'^((7, B, t) given 
in Theorem 14. II Then we find 

{m^\^{xM = (a/3i)-il|||,||<,,}a; 

r„i„v] . ^ «/35 . /IkiP \ a(/32 + fc/33) ^ (5-21) 

[7^ J3 5'(x,0o) - -^^^{\\A\<r} \^ -kj (1{||,||<,} -7) . 

Together with (j5.19p . this immediately yields the expansion for ^/njln and the expansion for 
[Pn - r/a) with 



Ai 
A2 



(35 



2a7/36 2fc7Q;'^ ' 
a(/32 + fc/33 + fc/35) r3 



27/36 2fc7a3 /Sg ' 



a(/32 + fc/33) 7 , /, 2 2\ 
2^-^+2fc^^'" 
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To obtain the expansion for the covariance estimator, note that P satisfies the conditions of 
Theorem 4.2 in [3]. This means r„ — > al with probabihty one, so that 

E„ - a'l = (f „ + a/)(r„ - al) = 2a(r„ - al) + o(l), 

with probabihty one. Hence, from (j5.18p we obtain 

V^(s„-a2/) =__^^[z?»-]^xi,(X„0o) + op(l), (5.22) 



where 



''^^'^'-^'^'\,.,<_rA^^-k].I (5.23) 



27/?2/3( 



'6 



a (3 A 



+ ^(l{INl<r}-7)-/. 



This yields the expansion for ^/n ( I]„ — a^/ ) with 



2a ka^{f33(3G - (3^(35) a^p. 



«2 



K3 



K4 



/32 7/^2/36 7/36 ^7' 

PiP^ - PsPe ^ aP2 + 27 

7/32/36 A;7q;/32 
2 

a/32' 



Proof of Theorem 14. 2t The expansion for y/nfin given in CoroUary 14. 1[ together with the 
fact that E{||Xi|| < r}Xi = 0, yields that ^/n'jj.n is asymptotically normal with mean zero and 
covariance matrix ^ 

^2e{||Xi|| < r}XiX[ = yE{||Xi|| < r}\\X,f ■ I. 
Since tt = — (q;/3i)~^, together with (14. 3|) . we find 

r=^E{\\X,\\<r}\\X,r= 



which proves part (ii). To prove (iii), first note that from Corollary [131 it follows that 

1 . /XX' \ 



where £{y) = K3l{||j,||<r}?; and m{y) — l{||j,||<r}(Ki + K2J/ ) + K4. Note that according to (|4.3 

E£(||Xi||) = -^E{||Xi|| < r}\\X^r = 

a/32 fh 
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so that E[^(||Xi||) + fcm(||Xi|l)] = 0. Since also E^2(||Xi|l) < oo and ETn^i\\Xi\\) < oo, it follows 
from Lemma 5 in |TU] , that the sum on the right hand side of (j5.24[) is asymptotically normal with 
mean zero and covariance matrix cri(/ + Ck,k) + cr2vec(/)vec(/)', where 

'^^ k{k + 2) ' 

^2 = ^^^'™^ +IEm^(ll^i||) + ^E£(||Xi||M||Xi||). 
If we fill in the expressions for ^(||Xi|| and rn(||Xi||), we get 

+ ^2(ki + K4)k2 + + K4)^ El{;||Xj|[<r}||Xi|p + Ki(ki + 2k4)E1{ || <r} + K4. 

Substituting the expressions for ki,K2,K4 given in Corollary 14.31 together with (|4.3p and (|4.4|) 
proves (iii). For part (iv) note that 

E [Ail{||x,||<r}ll^»lP + A2l{||x,||<r} + A3] = Xik-fa^ + X2I + A3 = 0. 

Therefore, from the expansion given in Corollarv l4.11 it follows that ^Jn{pn — r/a) is asymptotically 
normal with variance 

= E (Ail{||jc^||<r}||Xi|p + \2t{\\x,\\<r] + A3)^ 

= A2El{||x,||<r}||^^f + Al(A2 + A3)El{||x,||<.}||X,||2 + A2(A2 + A3)El{||x,||<.} + \l 



Substituting the expressions for A2, A3 given in Corollary 14. 31 together with ()4.3|) and ()4.4p proves 
(iv). Finally, for part (i), first note that according to Theorem 5.1 in [3], I]„ and p„ are 
mutually asymptotically normal. Hence, it suffices to prove that the quantities considered in part 
(i) are asymptotically uncorrelated. However, this follows directly from the expansions given in 
CoroUarv 14. II together with the symmetry properties of spherically symmetric densities. I 

Proof of Corollary 14. 3t According to Theorem 1 in ,2 , the MCD functional 6*0 = (/io,ro,po) 
as defined in (12. 6p is unique, and since P has a density, all conditions of Theorem 5.2 in [3] 
are satisfied. It follows from this theorem that the influence function for the functional Q{P) = 
{^i{P),T{P),p{P)), where T{Pf = is given by 

IF(x; e, P) = -A'iOoy'^ix, 0o), (5.25) 

where ^ is defined in (|2.7p . The expressions for IF(a;; /x, P) and IF (a;; p, P) follow directly from (|5.2ip . 
To obtain the influence function for the covariance functional, first note that according to the con- 
tinuity of the MCD functional, r(Pe,^) ^ r(P) = al, as e i 0, where P^.x = (1 - £)P + e4- This 
means that 

i](p,,,) - s(P) = (r(p,,,) + r(p)) (r(p,,, - r(p)) = 2a(r(p,,,) - r(p)) + o(e). 

It follows that 

IF(a;; E, P) = 2a • IF(a;; F, P) = -2a [P>'"^] ^ *(a;, 0o)- 
The expression then follows from (|5.23p . ■ 
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